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ABSTRACT 

We study scaling relations of compressible strongly magnetized turbulence. We find a good corre- 
spondence of our results with the Fleck (1996) model of compressible hydrodynamic turbulence. In 
particular, we find that the density-weighted velocity, i.e. u = p^^'^v, proposed in Kritsuk et al. (2007) 
obeys the Kolmogorov scaling, i.e. Eu{k) ^ k^^/^ for the high Mach number turbulence. Similarly, we 
find that the exponents of the third order structure functions for u stay equal to unity for the all the 
Mach numbers studied. The scaling of higher order correlations obeys the She-Leveque (1994) scalings 
corresponding to the two-dimensional dissipative structures, and this result does not change with the 
Mach number either. In contrast to v which exhibits different scaling parallel and perpendicular to 
the local magnetic field, the scaling of u is similar in both directions. In addition, we find that the 
peaks of density create a hierarchy in which both physical and column densities decrease with the 
scale in accordance to the Fleck (1996) predictions. This hierarchy can be related ubiquitous small 
ionized and neutral structures (SINS) in the interstellar gas. We believe that studies of statistics of 
the column density peaks can provide both consistency check for the turbulence velocity studies and 
insight into supersonic turbulence, when the velocity information is not available. 
Subject headings: ISM: structure — MHD — turbulence 



1. INTRODUCTION 

The interstellar medium (ISM) is a highly compressible 
turbulent, magnetized fluid, exhibiting density fluctua- 
tions on all observable scales. It has been long realized by 
many researchers that incompressible hydrodynamic, i.e. 
Kol mogorov, description is ina, dequate for such a medium 
(see lElmegreen fc Scalo |[2004 for review). Scaling rela- 
tions, if they were obtained for the interstellar gas, would 
be very helpful for addressing many problems, including 
the evolution of molecular clouds and star formation. 

Attempts to include effects of compressibility into the 
interstellar turbulenc e description can be d ated as far 
back as the work by Ivon Weizsackeij ()1951[ ). There a 
simple model based on a hierarchy of clouds was pre- 
sented. According to this picture every large cloud con- 
sists a certain number of smaller clouds , which contain 
even s maller clouds. For such a model Ivon Weizsackeij 
()1951h proposed a relation between subsequent levels of 
hierarchy 
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(1) 



where p^j is the average density inside a cloud at level 
i/, li, is the mean size of that cloud, 3 is the number of 
dimensions, and a is constant that reflects the degree of 
compression at each level v. 

The Kolmogorov energy spectrum (~ fc~^/^) follows 
from the assumption of a constant spe cific energy trans- 
fer rate e ~ jiljv). iLighthill! (jl955f ) pointed out that, 
in a compressible fluid, the volume energy transfer rate 
is constant in a statistical steady state 



sv = pe ^ pv'^ /(l/v) = pv^ jl. 



(2) 



In an imp ortant, but not sufficiently appreciated work, 
iFleckl (|1996il (henceforth, F96) incorporated above hier- 
archical model with energy transfer in compressible fluid 
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to obtain the scaling relations for compressible turbu- 
lence. By combining the equations ([1]) and ([2]) iFleckl 
(|1996f ) presented the following set of scaling relations in 
terms of the degree of compression a: 

where Ni and Mi are, respectively, the column density 
of the fluctuation with the scale I and the mass of the 
cloud of size I. The fluctuations of velocities in F96 model 
entail the spectrum of velocities £'(fc) ~ fc~^ / ^~^". 

In the spirit of F96 model, .Kritsuk et al.l ()2007D pro- 
posed to use the density-weighted velocity u = p^/^v as a 
new quantity, for which the Kolmogorov scaling for sec- 
ond order structure functions (SFs) can be restored (see 
Eq. [2]) in compressible hydrodynamic turbulence. Their 
hydrodynamic simulations provided for u spectrum close 
to Kolmogorov and the third order structure function 
that scales in proportion to distance. 

Will the F96 model be valid for compressible strongly 
magnetized turbulence? This is the major question that 
we address in this paper. 

2. NUMERICAL MODELING 

We used an second-order-accu rate essentially nonoscil- 
latory (ENO) scheme (see ICho fc Lazari an 2002J; 
Kpwal et al. 2002, for details) to solve the ideal isother- 
mal magnetohydrodynamic (MHD) equations in a pe- 
riodic box with maintaining the V • B = constraint 
numerically. We drove the turbulence at wave scale 
k ~ 2.5 (2.5 times smaller than the size of the box) 
using a random solenoidal large-scale driving accelera- 
tion. This scale defines the injection scale in our models. 
The rms velocity 5v is maintained to be approximately 
unity, so that v can be viewed as the velocity measured in 

1/2 

units of the rms velocity of the system and B / (47rp) ' 
as the Alfven velocity in the same units. The time t is 
in units of the large eddy turnover time (^ L/6v) and 
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Fig. 1. — Spectra of velocity and density- weighted velocity 
(dashed and solid lines, respectively) for super- and subsonic mod- 
els (big and small plots, respectively). Spectra are compensated 
by fcS/a. 

the length in units of L, the scale of the energy injec- 
tion. The magnetic field consists of the uniform back- 
ground field and a fluctuating field: B = Bcxt + b. Ini- 
tially 6 = 0. We use units in which the Alfven speed 

1/2 

VA = Bext/ (47rp) ' — 1 and p — 1 initially. The values 
of Bcxt have been chosen to be similar to those observed 
in the ISM turbulence. For our calculations we assumed 
that Scxt/ (47rp)^''^ ~ SB/ {Anp)^^'^ ^ Sv. In this case, 
the sound speed is the controlling parameter, and ba- 
sically two regimes can exist: supersonic and subsonic. 
Note that within our model, supersonic means low /?, i.e. 
the magnetic pressure dominates, and subsonic means 
high P, i.e. the gas pressure dominates. 

We present results for selected 3D numerical experi- 
ments of compressible MHD turbulence with a strong 
magnetic field for sonic Mach numbers Ais between « 0.7 
and 7. The Alfvenic Mach number Ma ^ 0.7. Mach 
numbers are defined as the mean value of the ratio of the 
absolute value of the local velocity v to the local value 
of the characteristic speed Cs or va (for the sonic and 
Alfvenic Mach number, respectively). To study effects of 
magnetization we also performed superAlfvenic experi- 
ments with A4a ^ "2- AH models were calculated with 
the resolution 512"^ up to 6 dynamical times. 

3. RESULTS 
3.1. Kolmogorov Scalings for Supersonic Flows 

In Figure [1] we present the spectra for velocity v and 
density-weighted velocity u = p^^^v for two subAlfvenic 
highly magnetized, i.e. low-/3, models. Naturally, for 
subsonic model the differences between spectra for v 
and u are marginal and both spectra correspond to Kol- 
mogorov's fc~^/'^ scaling (see subplot in Fig. [1]). How- 
ever, we can see that for the supersonic case, the veloc- 
ity spectrum gets steeper. The steepening corresponds 
to a « 0.23 (from Ey ~ k'^^^-^"). At the same time, 
the spectrum of it matches well the Kolmogorov slope. 

In the original Kolmogorov theory ()Kolmogorovlll94ll 
hereafter K41) it was shown that the spectral index of the 
third order structure function (SF), e.g. structure func- 
tion of u, S^u\l) = {\u {r + l)-u (r) \^) - should 
be equal 1, i.e. ^3 = 1- Figure O we show the SFs 
of the third order for velocity and for density-weighted 
velocity for supersonic model. We checked, that for sub- 
sonic motions, for both v and u, the index ^3 is indeed, 
close to unity. For the supersonic case, (.3 increases with 
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Fig. 2. — SFs of the third order for v (dashed line) and u (solid 
line) compensated hy for two MHD turbulence models: sub- 
sonic and supersonic (top and bottom plots, respectively). Dotted 
lines show the Kolmogorov scaling S3 (J) ~ I. Two dotted vertical 
lines bound the intertial range. 

the Mach number for v, but stays the unity for u (see 
Fig. [2]). This suggests that the Kolmogorov universality 
is preserved for supersonic MHD turbulence when den- 
sity weighting is applied. 

3.2. She-Leveque Intermittency Model 

A proper description of turbulence requires higher mo- 
ments (see Lazarian 2006a, for review). Those charac- 
terize intermittency, which is in the original K41 model 
is not accounted for. A substantial progress in under- 
standing turbulence intermittenc y is related to a dis- 
covery bv IShe fc Leveauel (|1994 hereafter SL94), who 
found a simple form for the scaling of the spectral in- 
dex of higher order longitudinal correlations S'^P-'(Z) = 
(I [u{r + l) -u{r)] ■ i\P) - l^". While in K41 model 
(p = p/3, SL94 provides 
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(4) 



where g is related to the scaling of the velocity vi ~ Z^^^ , 
X is related to the energy cascade rate ~ l~^ and D 
is the dimension of the dissipative structures. In hydro- 
dynamic incompressible turbulence, we have 5 = 3 and 
X = 2/3. For MHD turbulence the dissipation happens 
in current sheets, which are two-dimensi onal dissipative 
struc tures, corresponding to D — 2 (Miiller fc Biskanrpl 
I2OOOI) . Thus, for subsonic MHD turbulence we expect 
Cp = I + 1 - (l/3)P/3 for both velocity and the den- 
sity weighted velocity. This is what we actually observe 
in Figure [3] (see subpanel). The same scaling, however, 
is preserved for u for supersonic magnetized turbulence, 
which indicate, that, unlike, velocities, u exhibits uni- 
versal intermittency. 

3.3. Anisotropics Induced by Magnetic Field 

Magnetic field is known to induc e aniso tropics 
of compressible MHD turbulence (see iHigdonI |l984). 
Anisotropy increasing with the de crease of scale was 
predicted for Alfvenic motions by iGoldreich fc Sidhail 
(1995:, henceforth GS95, see also iLithwick fc G oldrcichl. 
l200li) and confirme d numerically for compressible MHD 
in lCho fc Lazarian! (2002, 2003). 

For supersonic motions Figure 0] shows that the scal- 
ings for V are much steeper in both directions than 
those predicted by GS95 model (1.215 ± 0.007 and 
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Fig. 3. — Scaling exponents for v and u for subsonic (subplot) 
and supersonic models. The plots present unnormalizcd values 
of the scaling exponents obtained directly from SFs by fitting the 
relation S^^\l) = al^v within t he intertial range , i.e. without using 
the extended self-similarity (Ben zi et al. IIT993I ). 
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Fig. 4. — SFs of the second order for v and u in the local ref- 
erence frame for the supersonic experiment. SFs for velocity scale 
as ~ Zl-215±0.007 and ~ z0.S74±0.004 foj, parallel and perpendicu- 
lar directions, respectively. SFs for u scale as ~ ;0.882±0.006 ^nd 
~ ;0.744±0.003 £qj. II and _L directions, respectively. 

0.874 ± 0.004 for || and ± directions, respectively). How- 
ever, those slopes still give a close to GS95 anisotropy 
(/|| ~ ^o.7i8±o.oo2^^ which is indicative of the dominance 

of the Alfvenic ("incompressible") motions. Note, that 
in Figure 2] the SFs are obtained in the system of ref- 
erence of the local magnetic field, i.e. the field on the 



scales of the fluctuations under study. S^p and S^^' de- 
note second order SFs parallel and perpendicular to local 
magnetic field, respectively. 

For u the slopes are significantly smaller (0.882 ±0.006 
and 0.744 ± 0.003 for || and ± directions, respectively). 
The SF in perpendicular direction scales more like incom- 
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pressible motions, i.e. S^l'' ^ Z^''^. The slope of 5|| 
u is smaller than the corresponding one for v resulting in 
the reduced degree of anisotropy (/|| ^ /0-843±o.oo4-j^ Intu- 
itively, this can be understood in terms of dense clamps 
strongly distorting magnetic field as they move in respect 
to magnetized fluid. 

3.4. Statistics of Column Density Peaks 

Our results for velocity show that our simulations of 
strongly magnetized turbulence provide a ~ 0.23 for 
Ais ~ 7. According to F96 model (see Eq. [3|) this sug- 
gests the existence of a rising spectrum of density fluctu- 
ations within the hierarchy of density clumps. We try to 
make our study more related to observations which usu- 
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Fig. 5. — Scaling relations for the column density A' for three 
models of turbulence: subsonic {A4s ~ 0.7) and supersonic {A4s ^ 
2 and 7, see legend). 

ally measure densities integrated along the line of sight, 
i.e. column densities, or alternatively study the hierarchy 
of observed clump masses (see Eq. [3]) . 

F96 model assumes the existence of an infinitely ex- 
tended hierarchy. In our computations the structures are 
generated by turbulence at scales less than the scale of 
the computational box. Therefore F96 scaling relations 
(in Eq. [3|) should be modified as follows 
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Our procedure of obtaining the scaling relation from 
column densi t y ma ps can is similar to that, e.g. in 
iKritsuk et al.l ()2007[ ). with the difference that they dealt 
with 3D data, while we deal with 2D data. First, we seek 
for a local maximum of column density. Then we calcu- 
late the average column density within concentric boxes 
with gradually increasing the size I. In case of determin- 
ing the relation for Mi , instead of averaging we apply the 
integration over the boxes. Naturally, the results should 
correspond to each other. 

In Figure [5] we present an example of scaling relation 
for column density for three models of turbulence with 
AAs ^ 0.7, 2, and 7. One can note, that the relation be- 
comes more steep with the sonic Mach number within the 
intertial range. The fractal dimensi ons can be calculate d 
from the relation Dm = 3-1-7 (see IKritsuk et al.l[200l . 
where 7 = —3a is a slope estimated from the plot within 
the intertial range. For our models the fractal dimension 
ranges from Dm ~ 2.5 for the highly supersonic models 
to Dm — 2.9 for subsonic model. Respectively, the com- 
pressibility coefficients for presented models a ~ 0.04 for 
Ms ~ 0.7 to a ~ 0.19 for Ms 7. The latter roughly 
consistent with a obtained for the velocity SF measure- 
ment in i )3.1l The differences are probably due to insuf- 
ficient statistics of rather rare high peaks. In general, 
the filling factor of peak decreases with the height of the 
peak. 

3.5. Variations of Scalings Induced by Fluid 
Magnetization 

What is the effect of magnetic field on the u-scaling? 
The spectra, third and higher moments of correlations 
obtained for our superAlfvenic simulations with Ma 2 
happen to be very similar to t he case of strongly m ag- 
nctized turbulence. Results in IKritsuk et al.l (|2007f ) on 
the Kolmogorov spectrum of u obtained for pure hydro 
(i.e. Ma ^ 00). Our results indicate that, unlike veloc- 
ity, u is much less affected by magnetic field. Naturally, 
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in super- Alfvenic turbulence the anisotropies induced by 
magnetic field are not observed at large scales within the 
intertial range (cf. the last paragraph of §31) • 

4. ASTROPHYSICAL IMPLICATIONS 

Dependence of a on the extend of inertial range — If we 
combine several facts together, namely, (a) that a is 
a function of Mach number, (b) that the maximum of 
density correspond to the dissipation scale, e.g. shock 
thickness scale Idissj (c) that the amplitude of density 
in peaks scales as the mean density times A^^, we have 
to conclude that as the inertial range from the injection 
scale linj to Idiss increases, for a given Mach number, 
a should decrease. Connecting these facts we get the 
following relations, ppeak ^ Al? ^ {Imj /Idiss)^" , which 
gives the dependence of a on Unj /Idiss and Ms, namely, 
a logAis/ ^og{li„j /Idiss)- An interesting consequence 
of this would be a prediction of Kolniogorov scaling for 
supersonic velocities when the injection and dissipation 
ranges are infinitely separat ed. Consequent l y, the steeper 
velocity spectra reported in lPadoan et al. I ( 20071 ) can be 
interpreted as an indication of a limited inertial range. 
Further research justifying such conclusions is required, 
however. 

SINS of supersonic turbulence — Ubiquitous small ion- 
ized and neutral structure s (SINS) a re observed in in- 
terstellar medium (see H eiles. I Il997l ). Their nature is 
extremely puzzling if one thinks in terms of Kolmogorov 
scalings for density fluctuations. The fact that the spec- 
trum of fluctuations of density in supersonic turbulence 
is shallower that the Kolmogorov one is well-known (see 
[Kowal ct al. 2007, and references therein). However, just 
the difference in slope cannot explain the really dramatic 
variations in column densities observed. The present 
paper provides a different outlook at the problem of 
SINS. Wc sec that, while low amplitude density fluctua- 
tions ex hibit Kolmogorov scaling (Bcrcsnyak ct al. 200^ 
iKowal et al.l 120071 ) . high peaks of density correspond to 
a rising spectrum of fluctuations. Thus, observing super- 
sonic turbulence at small scales, we shall most frequently 
observe small amplitude fluctuations corresponding to 
Kolmogorov-like spectrum of density fluctuations. Occa- 
sionally, but inevitably, one will encounter isolated high 
density peaks. An alternative mechanism for getting in- 



frequent lar ge density fluctua tions over small scales is 
presented in iLazarian I (|2006bD and is related to current 
sheets in viscosity-damped regime of MHD turbulence. 

Clumps and star formation — Interstellar medium is 
known to be clumpy. Frequently clumps in molecular 
clouds are associated with the action of gravity. Our 
study shows that supersonic turbulence tend to produce 
small very dense clumps. If such clumps happen to get 
Jean's mass, they can form stars. Therefore, star forma- 
tion is inevitable in supersonic turbulence. However, the 
efficiency of star formation is expected to be low, as the 
filling factor of peaks decreases with the increase of the 
peak height. Inhibiting of star formation via shearing 
may dominate in terms of influencing of star-formation 
efficiencies. 

5. SUMMARY 

In paper above we have studied the scaling of super- 
sonic MHD turbulence. We found that: 

• Fleck 1996 model is applicable to strongly magne- 
tized compressible turbulence. 

• Spectra and structure functions of density- 
weighted velocities are consistent with Kolmogorov 
expectations. 

• Intermittency of density-weighted velocity can be 
well described by the She-Leveque model with the 
dimension of dissipative structures equal 2. 

• Strongly magnetized supersonic turbulence demon- 
strate lower degree of anisotropy if described using 
the density-weighted velocity. 

• The high peaks of column density exhibit increase 
of the density with the decrease of scale, which may 
be relevant to the explanation of SINS. 
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